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ABSTRACT 

This study aimed to derived an equation to estimate the monthly rainfall 
for Catarman, Northern Samar.The observed monthly rainfall data for 
Catarman N. Samar, Catbalogan Samar, Legazpi City and Masbate were 
obtained from the Philippine Atmospheric Geographical Astronomical 
Services Administration (PAGASA). The monthly rainfall records of the 
three (3) neighboring stations (Catbalogan, Legazpi, Masbate) were used to 
identify which of the existing rainfall prediction methods, namely, Normal 
Ratio Method, Distance Power Method and Multi Linear Regression Method 
is the basis in the development of a new equation. 


The accuracy by which the existing methods predict the observed monthly 
rainfall in Catarman was evaluated using T-test for correlated samples and 
the Pearson’s Correlation Coefficient. Since none of the methods produced 
estimates nearest to the observed monthly rainfall in Catarman, an 
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equation has been derived: 
mee 13 NRA cn Pidi 7 
T i 


Using the neighboring stations of Catbalogan, Legazpi, and Masbate the 


equation became: 


y= 21 { L 13 NRa DcCATBALOGANP CATBALOGAN 4 DLEGAZPIPLEGAZPI 4 
25 20 Dr NRCATBALOGAN NRLEGAZPI 
DMASBATEP MASBATE \" 
NRMASBATE 


Keywords: Monthly Rainfall, PAGASA, Normal ratio method, distance power 


method, Multi linear regression method 


INTRODUCTION 

Rainfall data are one of the most important and frequently 
used hydrological data in water resources planning and 
other climatological analyses. But rainfall is affected by 
climate change which is being experienced nowadays. One 
of its tremendous effects is the variably unexpected 
rainfall intensity. This phenomenon cannot be explained 
by natural variation under normal conditions. El Nino and 
La Niffa phenomena affect the rainfall pattern and 
intensity. 


Pajuelas (2000) observed that rainfall distribution over 
Space and time during either El Nino or La Nina event is 
highly variable and causes major economic losses due to 
disastrous typhoons, floods, and droughts. It is during 
these extreme rainfall events that rainfall characteristics 
such as rainfall depth and rainfall intensity should be 
monitored and measured so that such data may be used 
for weather forecasting and planning purposes. 


Rainfall data is very important and useful in the planning 
of agricultural development, water resources and disaster 
mitigation projects and programs. Hence, functional rain 


URL: 
www.ijtsrd.com/papers /ijtsrd35875.pdf 


Copyright © 2020 by author(s) and 
International Journal of Trend _ in 


Scientific Research @) 
and Development AGA - 
Journal. This is an 7 : 

Open Access article distributed under 
the terms of the Creative Commons 
Attribution License (CC BY 4.0) 
(http: //creativecommons.org/licenses/by/4.0) 





gages should be installed in all municipalities of the 


country or much better in every barangay. The more rain 
gages there are, the more accurate is the measurement of 
rainfall as well as the estimation of the average rainfall 
over an area where a development project is being 
proposed. 


The problem, however, on rainfall monitoring and 
measurement is that, first and foremost, there is one rain 
gage for every 30,000 hectares (Linsley, et al., 1982). In 
the Philippines, a major problem is the continuous 
functionality of rain gages because many of them are 
destroyed during the typhoon months (September - 
December) and they are not immediately repaired or 
replaced by a functional one. 


The continuity of precipitation data over time is also very 
important. For the rainfall forecasting to be credible, this 
should be based on a 30-year rainfall record. Almost all 
rainfall data have gaps not only due to loss of records 
because of fire accidents, wars, typhoons, floods, etc. but 
also due to occasional interruption abnormal function of 
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the instruments used. Hence, numerous techniques were 
developed to estimate the missing precipitation data. 
Pajuelas (2000) stated that estimation of missing data is 
known as the first stage of most climatological, 
environmental and _ hydrological studies. Spatial 
interpolation of precipitation data for imputation of 
missing records is an essential and crucial step in the 
development of continuous precipitation data without any 
gaps needed for hydrologic analysis, modeling and design 
(Teegavarapu, 2011). 


Developing countries, like the Philippines, are often faced 
with this problem. There are many methods for estimation 
of missing rainfall data. According to Hasan, et 
al.(2013)spatial interpolation techniques are widely used 
methods for filling the gaps in daily rainfall series through 
estimating the unknown rainfall amount for a point from 
the known data of adjacent stations as proven by the 
studies of previous researchers. Paulhus and Kohler 
(1952) explored two methods of interpolation, the 
normal-ratio and 3-station-average, to fill the missing 
values in monthly rainfall data. The Inverse Distance 
Weighting (IDW) methods estimate the rainfall amount of 
a location as a weighted average of the rainfall amount of 
adjacent stations and the weights are considered as a 
function of the distances (Teegavarapu and Chandramouli, 
2005). The correlation coefficients between data Series are 
also explored to estimate the weights (Ahrens, 2006). 
Regression based methods are also used for estimating 
missing precipitation values (Lo Presti et. al., 2010). 
Regression models consider climate data, elevation, 
topography, proximity to coastal area etc. as explanatory 
variables to estimate missing rainfall series of a station 
(Daly et. al., 1994). 


Most of these methods in estimating missing rainfall data, 
however, were developed applying rainfall data from 
other countries, especially the USA, United Kingdom, 
Australia and other highly developed ones. The rainfall 
pattern, frequency and intensity in the Philippines is 
different from those countries. Moreover, the province of 
Northern Samar has different rainfall characteristics 
compared to the other provinces of the Philippines and so 
Catarman has different rainfall characteristics from the 
other towns of the province of Northern Samar. 


These problematic realities challenged the researcher to 
develop an equation on estimating the monthly rainfall for 
Catarman, Northern Samar. 


This study aimed to develop an equation to estimate a 
missing monthly rainfall specifically for Catarman, 
Northern Samar. 


METHODOLOGY 

This study was conducted in the Municipality of Catarman, 
the capital town of the Province of Northern Samar. It is a 
first-class municipality and is the largest town in terms of 
land area and population in the province. It lies on the 
northern part of the Samar Island. It is bounded to the East 
by Mondragon, to the west by Bobon, to the south by Lope 
de Vega and to the North by the Philippine Sea. 


Catarman has a tropical climate. There is significant 
rainfall throughout the year in Northern Samar. Even the 


driest month still has lot of rainfall (Pajuelas, 2000). The 
Philippine Atmospheric, Geophysical and Astronomical 
Services Administration (PAGASA) Regional Office No. VIII 
has its Satellite Office in Catarman, Northern Samar 
situated at Brgy. Dalakit of this town. 


This study employed the quantitative correlational study, 
a descriptive research design. In this study, the correlation 
of the rainfall data of Catarman and of the six identified 
neighboring stations, the Catbalogan, Tacloban, Borongan, 
Sorsogon, Legazpi and Masbate synoptic stations, was 
done to determine whether the data of the latter can be 
used to estimate for the missing monthly rainfall of the 
former. 


As regards data gathering, the observed daily rainfall data 
for the period January 1, 1986 to December 31, 2015, at 
the Catarman Synoptic Station, Catarman, Northern Samar 
were personally obtained by the researcher from Mr. 
Felion C. Corona, the Chief Meteorological Officer of the 
Said station. 


For the six neighboring stations which were originally 
included in the study, the observed daily rainfall data for 
the same period were obtained from the Climatology and 
Agro meteorology Division ofthePhilippine Atmospheric 
Geophysical and Astronomical Services Administration 
(PAGASA) Central Office, Diliman, Quezon City. 


The researcher looked for methods of estimating missing 
monthly rainfall in the different books and journal articles 
reviewed. Based on the data requirements and level of 
accuracy, the researcher has chosen three methods such 
as the Normal Ratio Method, Distance Power Method and 
Multi Linear Regression Method. The Normal Ratio 
Method can be used only if the available daily rainfall data 
for each of the stations involved in the study runs for 
thirty-years, from January 1, 1986 to December 31, 2015. 
For the Distance Power Method and Multi Linear 
Regression, the available rainfall data should be at least for 
ten (10) continuous years. 


The stations which qualified for the Normal Ratio Method 
are the Catbalogan, Legazpi and Masbate stations while 
those qualified for Distance Power and Multi Linear 
Regression Methods are Catbalogan, Legazpi, Masbate, 
Tacloban and Borongan Stations. For purposes of 
correlation and consistency of the data, the researcher 
decided to consider only the stations having the complete 
thirty years rainfall data set. Thus, only the stations of 
Catbalogan, Legaspi and Masbate were used in this study. 


The data obtained were examined as to _ their 
completeness and sufficiency for the application of the 
three identified methods to be used for the estimation of 
missing rainfall. Based on this evaluation, only three 
neighboring stations qualified and these are _ the 
Catbalogan, Legazpi and Masbate Synoptic Stations. 


This study utilized the daily rainfall from the Philippine 
Atmospheric Geophysical and Astronomical Services 
Administration (PAGASA) of Catarman and the three 
neighboring synoptic rain gage stations, namely, 
Catbalogan, Legazpi and Masbate, instead of the original 
plan to use six stations because the Sorsogon, Borongan 
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and Tacloban stations were excluded by reason of 
incomplete data series. 

The daily rainfall gathered from the Catarman, Catbalogan, 
Legazpi and Masbate Synoptic Stations were summed up 
to find the monthly rainfall. 


A Double Mass Analysis was performed to determine the 
homogeneity of the rainfall data set from Catarman Station 
and the rainfall data sets from the identified three 
neighboring stations. A break in the slope of the graph 
indicated a change in the precipitation regime in a 
particular station and the record of annual rainfall was 
adjusted using the formula 


Pa= -< Po (equation 1) 


where: 

Pa = adjusted precipitation 

Po = observed precipitation 

ba = slope of graph to which records are adjusted 
bo = slope of graph at time Po was observed 


As shown in Figure 1, the observed rainfall of the three 
neighboring stations, namely, Catbalogan, Legazpi and 
Masbate and that of the test station (Catarman) possess 
homogeneity and consistency. 
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Figure 1 Double mass analysis curve for the four (4) 
sites 


Having proven that the rainfall data for Catarman station 
and the three neighboring stations are homogeneous and 
consistent, the researcher started the computation of the 
estimated monthly rainfall for Catarman station using the 
monthly rainfall data from the three neighboring stations 
by the following formulas: 


dt copy NR . 
A. Normal Ratio Method P, = et ao (equation 2) 

L 
where 
Px = the missing rainfall for Catarman station at a 
particular month 


P; = the rainfall for the same month at any of the 
neighboring station either Catbalogan, Legazpi or Masbate 
NRa = the normal monthly rainfall value for Catarman 
station 


NR; = the normal monthly rainfall value for any of the 
neighboring station either Catbalogan, Legazpi or Masbate 


NR = normal monthly (mean of thirty 30 years of monthly 
rainfall data) 


n = number of surrounding stations whose data are used 
in the estimation 


B. Distance Power Method 


Mbase Pi,j 
viz 


“i (equation-3) 


Py = yMbase 1 
i 


where: 
Px = missing monthly rainfall at Catarman station at time 


assy) 


J 


Pij = observed monthly rainfall at any of the neighboring 


as) 


Stations at time °j 


Dj = distance of Catarman station from any of the 
neighboring stations 


Mbase = number of neighboring stations taken into account; 
and 


b = 2, power of distance D used for weighting rainfall 
values at individual station. 


and the distance Dj is computed by the formula: 
Di = (x- xi? + (y- yi)? + (2-2) 


where: 

xX, y, Z = are the coordinates of the Catarman Station; and 

Xi, yi, Zi = are the coordinates of any of the neighboring 
Stations. 


3. Multiple Linear Regression Method 
y = bo + bi xi + b2 x2 + ... + DnXn 


where: 

y = missing monthly rainfall of Catarman station 

X1, X2, .... Xn= monthly rainfall at the neighboring stations; 
and 

bo, bi, bz, .... bn = partial regression coefficients determined 
by the statistical package using the monthly rainfall data 
from the neighboring stations. 


The statistical procedures employed in this study were the 

following: 

1. “T” test for Correlated Samples (to determine if there is 
significant difference between the estimated missing 
monthly rainfall (at Catarman station) using the three 
methods and the observed monthly rainfall at 
Catarman station) 


t= ce (equation-4) 


D 
¥ D2 oa OD) 
n(n—-1) 
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X = monthly rainfall in any neighboring station included in 


where: the study; 

D = main difference between observed monthly rainfall at Y = the observed monthly rainfall data at Catarman 
Catarman station and the estimated missing monthly Station; and 

rainfall by any of the three methods employed n = population size. 


n= sample size 
Finally, a new equation to estimate monthly rainfall data 


2. Pearson’s Correlation Coefficient (to determine the was developed. A method was selected as a trial equation 
degree of correlation between the observed and and was modified by incorporating new parameters of the 
estimated monthly rainfall for Catarman station) Stations (distance, elevation, etc.). The monthly rainfall for 

Catarman using the first trial equation is estimated and 
yxy -20GY) compared with the observed monthly rainfall in the 

R= — (equation-5) Catarman station by T test values and Pearson’s 

(zx2-2\( pve - ur) correlation coefficient. The trial equation was revised by 


adding empirical coefficients derived by mathematical 
transformation. This was repeated until the “t” value is 
very, very small and the correlation coefficient is almost 1. 


where: 
R = Pearson’s Correlation Coefficient; 


RESULTSAND DISCUSSION 
The normal monthly rainfall for Catarman, Catbalogan, Legazpi and Masbate stations were computed and presented in 
Table 1. 
Table 1 the normal monthly rainfall (mm) in the Catarman, Catbalogan, Legazpi, and Masbate stations. 
Month Catarman Catbalogan Legazpi Masbate 





Using these normal rainfall values, the monthly rainfall in Catarman, N. Samar was estimated by the Normal Ratio Method 
as shown in Tables 2a and 2b. 


Table 2a the estimated monthly rainfall in Catarman, N. Samar by the Normal Ratio Method (January to June) 
Estimated Monthly Rainfall 
January February March April May | June 
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| 2014 | 5280 | 66.9 | 307.4 | 244.6 | 33.2 | 226.5 | 


Table 2b the estimated monthly rainfall in Catarman, N. Samar by the Normal Ratio Method (July to December) 
Estimated Monthly Rainfall 
July August September October November December 


Using the geographical coordinates of Catarman, Catbalogan, Legazpi and Masbate stations, the computed distances are: 
Catarman-Catbalogan, 84,556.6087 m; Catarman-Legazpi, 197,338.71 mand Catarman-Masbate, 88,417.961 m. With these 
distances, the estimated monthly rainfall for Catarman were computed by the Distance Power Method as shown in Tables 
3a and 3b. 





Ma Ve 


Table 3a the estimated monthly rainfall in Catarman, N. Samar by the Distance Power Method (January to June) 
eee Estimated Monthly Rainfall 
January February March April May June 


80.4 
1990 
1992 


1995 
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| 2007 |_ 309.7 | 50.6 | 50.9 | 42.9 | 180.1 | 99.0 | 


Table 3b The estimated monthly rainfall in Catarman, N. Samar by the Distance Power Method (July to December) 
Estimated Monthly Rainfall 
July August September October November December 





Ma Ve 
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The coefficients used in the Multi Linear Regression Method were computed using a statistical package and presented in 
Table 4. 
Table 4 the computed coefficients used in the computation of monthly rainfall in Catarman N. Samar by the Multi 
Linear Regression Method 
Month Catbalogan Legazpi Masbate 


Tables 5a and 5b shows the estimated monthly rainfall using theMultI Linear Regression Method. 





Table 5a the estimated monthly rainfall in Catarman, N. Samar by the Multi Linear Regression Method (January to 
June) 
Estimated Monthly Rainfall 
January February March April May 
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Table 5b The estimated monthly rainfall in Catarman, N. Samar by the Multi Linear Regression Method (July to 
December) 
Estimated Monthly Rainfall 
July August September October November December 


iV 
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| 2002 | 173.5 | 2304 | 96.0 | 1001 | 178.2 | 1705 | 
| 2005 | 120.6 | 99.0 | 170.7 | 1729 | 1235 | 550.5 __| 
(2015 | 86.3 | 996 |. 62.6.0) .892 | 145.5 | 217.2 | 


The results of the “T” test for Correlated Samples (to determine if there is significant difference between the estimated 
monthly rainfall at Catarman station) using the three methods and the observed monthly rainfall at Catarman station) and 
the results of the Pearson’s Correlation Coefficient computation (to determine the degree of correlation between the 
estimated monthly rainfall by the three methods and the observed monthly rainfall for Catarman station are shown in 
Table 6. 





Table 6 the result of the T test for Correlated Samples and the Pearson’s Correlation Coefficient performed 
between the estimated monthly rainfall of the three methods and the observed rainfall in Catarman, N. Samar 
Observed Monthly Rainfall in Catarman Station versus 
ie Estimated Monthly Rainfall computed by the 


Conducted Normal Ratio Distance Power Multi Linear ee 


Method Method Regression Method 


T Test for Correla-ted Samples 1.645 
Pearson’s Correlation Coeffi-cient | 0.84 | 76 | 8A 





After several statistical and data transformation procedures, an equation for estimating monthly rainfall atCatarman, 
sores Samar has been derived, to wit: 








25 13 NR P,Dj 
[2 y}7 = = Log Dp, Lis 1 WR; (equation-6) 
simplifying 
21 13 NR P;D; ; 
y= =| Log Dp, dls 1 WR, |’ (equation-7) 
where: 


y = the monthly rainfall in Catarman 

P; = the rainfall for the same month at any of the neighboring station either Catbalogan, Legaspi or Masbate 

Dj = distance between Catarman and the neighboring stations (Catbalogan, Legazpi and Masbate). 

Dr= is the total distances between Catbalogan and Catarman, Legazpi and Catarman and Masbate and Catarman 
NRa= is the normal monthly rainfall value for Catarman station 

NRi= is the normal monthly rainfall value for any station 

n = number of neighboring stations 


Specifically, considering Catbalogan, Legazpi and Masbate as neighboring stations, the derived equation to estimate the 
monthly rainfall became: 


at { 13 NRa DCATBALOGANP CATBALOGAN is DLEGAZPIPLEGAZPI 4 DMASBATEP MASBATE ne 


y== (equation-8) 
25 20 Dr NRCATBALOGAN NRLEGAZPI NRMASBATE 
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where: 

4 = is the missing rainfall amount of Catarman 

Pcatbaloga = is the rainfall for the same month at 
Catbalogan station 

Plegazpi = is the rainfall for the same month at Legazpi 
station 

PMasbate = is the rainfall for the same month at Masbate 
station 

NRcatarman = iS the normal annual precipitation value for 


Catarman station 
NRcatbalogan = iS the normal annual precipitation value for 
Catarman station 


NRtegazpi = is the normal annual precipitation value for 
Catarman station 

NRmasbate = is the normal annual precipitation value for 
Catarman station 

Dcatbalogan = IS distance between Catbalogan and 
Catarman 

DLegazpi = is distance between Legazpi and Catarman 

DMasbate = is distance between Masbate and Catarman 

DtToTAL = is the total distances between Catbalogan 


and Catarman, Legazpi and Catarman, and 
Masbate and Catarman 


The normal annual rainfall, based on a 30 years rainfall 
data (1986 - 2015) for Catarman (Northern Samar) 
ranged from 154.9 mm (May) to 645.1 mm (December). 
For the same period, the normal annual rainfall in 
Catbalogan ranged 127.2 mm (April) to 339.9 mm 
(December); in Legazpi, it ranged from 176.6 mm (April) 
to 583.5 mm (December); and in Masbate, it ranged from 
53.7 mm (April) to 267.7 mm (December). These data 
shows that the three (3) neighboring stations have almost 
the same rainfall pattern with Catarman. These further 
shows that if there is an El Nino or La Nifa in the area, this 
is felt in all these four stations at the same time. Hence, the 
choice of the three neighboring stations is justified. 


To ascertain if the monthly rainfall at Catarman station 
can be predicted by the monthly rainfall in any of the 
neighboring stations, the coefficients of correlation 
between the monthly observed rainfall atCatarman station 
and any of the three identified neighboring stations were 
also determined. For Catarman and Catbalogan, the 
correlation coefficient was 0.737; for Catarmanand 
Legazpi, 0.622; and for Catarman and Masbate, 0.698. 
However, these correlation coefficients were very low. 


Based on the available rainfall record in the three (3) 
neighboring stations and level of accuracy, three methods 
to estimate missing monthly rainfall such as the Normal 
Ratio Method, Distance Power Method and Multi Linear 
Regression Method were identified. 


CONCLUSIONS 

In this study, the monthly rainfall in Catarman for 30 years 
was estimated by the three (3) methods. The estimated 
monthly rainfall using each of the three (3) methods was 
compared with the observed monthly rainfall in Catarman. 
The estimated monthly rainfall by the Normal Ratio 
Method has the lowest “t” value of - 0.3 which means that 
these estimated values are nearest to the observed 
monthly rainfall. The estimated monthly rainfall using the 
Normal Ratio and Multi Linear Regression Methods had 


correlation coefficients of 0.84 while the estimated 
monthly rainfall using the Distance Power Method had a 
correlation coefficient of 0.76 only. 


With these findings, the Normal Ratio Method was 
selected together with the Distance Power Method to be 
developed into a new equation for Catarman. The 
researcher believes that distance has some bearing on the 
rainfall pattern. An equation to estimate the monthly 
rainfall in Catarman has been derived. 


The “T” test for correlated samples on the observed 
monthly rainfall in Catarman versus the estimated 
monthly rainfall using the derived equation revealed a “t” 
value of 0.021 (<<<1.645) which means that the estimate 
monthly rainfall data is not statistically different from the 
observed monthly rainfall data in Catarman. The Pearson 
correlation coefficient for the same set of data was 
computed as 0.869 which means that the two data sets is 
highly correlated. The researcher further computed the 
estimate of missing annual rainfall for Catarman using the 
newly developed equation and gave a 0.986 degree of 
correlation which means that said equation can accurately 
estimate the missing rainfall data for the subject location. 
It can be concluded that the derived equation is acceptable 
and can estimate monthly rainfall in Catarman, N. Samar. 
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